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Abstract 

A method using an expansion of the four-body Yakubovsky wave function com¬ 
ponents onto the basis of the Faddeev-equation solutions for the two-cluster sub- 
Hamiltonian eigenfunctions is proposed. This expansion reduces the Yakubovsky 
differential equations to a system of coupled-channel equations for functions de¬ 
pending on the relative coordinates between the subsystems of the two-cluster 
partitions. On the basis of the resulting equations the four-nucleon bound- 
state problem and the zero-energy n-^H scattering problem are solved on the 
relatively small computer. 


1 Introduction 


The Faddeev-Yakubovsky approach provides the correct way for a treatment 
the few-body problem ii- It is based on the decomposition of the wave func¬ 
tion into components. This decomposition leads in the momentum representa¬ 
tion to integral equations. The solutions of these equations are uniquely defined 
by inhomogeneous terms in contrast to the Lippmann-Schwinger equations for 
which this holds only in the two-body case. In the coordinate representation, 
the decomposition generates differential equations |p which are easier to solve 
numerically because they involve only interparticle potentials and have simple 
boundary conditions |@]-|^]. 

For the three-body case these equations are widely used for practical calcula¬ 
tions. Unfortunately, already for the four-body systems, the numerical solution 
is very difficult and requires the most powerful supercomputer facilities |1C, 111. 
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Recently, two new methods were proposed to solve the differential Faddeev- 
Yakubovsky equations (DYE) in configuration space 0! 0) 0] • Both of them 
reduce the computational difficulties by several orders of magnitude, so that the 
four-body equations can be solved numerically on relatively small computers. 
One of the methods elaborated for the bound-state problem consists in 
the systematic exploitation of the tensor structure of the matrices appearing 
after discretization of DYE with spline expansion and orthogonal collocation 
methods. This technique looks very appropriate to calculations of the four-body 
bound states producing more than five significant digits for the four-nucleon 
binding energy. 

The key ingredient of the other method |12| is the partial diagonalization 
of four-body DYE by making use of the basis of the eigenfunctions of Ead- 
deev equations corresponding to subsystems of two-cluster fragmentations of 
the four-body system. Such expansions referred to conventionally as the cou¬ 
pled reaction channel (CRC) decomposition reduce DYE to a system of 
coupled equations for the functions depending only on the relative coordinates 
between the subsystems of two-cluster partitions. For this reason we name this 
procedure the cluster reduction. 

The purpose of this paper is to give the systematic description of the lat¬ 
ter method and to demonstrate its efficiency for numerical solutions of both 
the four-nucleon bound-state and scattering-state problems. The bound-state 
problem of the system nnpp ('^He) and the zero-energy n-^H scattering problem 
have been chosen for the actual calculations. The MT I-III and S3 potential 
models have been used as the N — N interactions 1^, 1^, 1^, 1^] . We would like 
to point out already here that all the calculations were performed on a personal 
computer with Intel DX-50 processor and 16 Mb of RAM. 

The paper is organized as follows. In the second section we recall those 
portions of the papers i i which are required to formulate the four-body 
bound-state and scattering-state problems on the base of DYE. In Sec. 3 the 
procedure of the partial diagonalization of DYE is described. There we give 
the appropriate definition of basis functions and discuss their main proper¬ 
ties of completeness and biorthogonality. The resulting CRC equations for the 
functions depending on the distance between the subsystems of the two-cluster 
partitions are used in the Sec. 4 for numerical calculations of the four-nucleon 
bound-state (^He) and zero-energy n-^H scattering problems. The last sec¬ 
tion summarizes the paper. We have collected some technical details in the 
Appendix. 


2 DYE formalism 

The four-particle wave function 'k should be decomposed into components in 
one to one correspondence to all chains of partitions. The chains consist of 
two-cluster partitions 02 {&■§■■, {ijk)l or {ij){kl)) and three-cluster partitions 03 
{e.g., {ij)kl) obeying the relation 03 C 02. The latter means that the partition 
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as can be obtained from the partition 02 by splitting up one subsystem. It is 
easy to see that there exist 18 chains of partitions for the four-particle system. 
The Yakubovsky wave-function components can be defined by the formulas 

[@> |] 

^a3a2 (E)Va, E ME)Vbs'^, (1) 

(b3f^a.3)Ca2 

where 

Ra3(E) = {Ho + Va, - E)-\ Ro{E) = {Ho - E)-\ 

Here Ho is the kinetic-energy operator in the c.m. frame and V^g stands for the 
two-particle potential acting inside the two particle subsystem of a partition 03 
{e.g., Kg = Vij with 03 = {ij)kl). 

If the function T is the solution of the Schrodinger equation 

{Ho + Y.Va3-E)^ = Q 

as 

then the components ^'0302 obey the Yakubovsky equations ||3|, 

{Ho + 1^3 ~ -^)'^^a 3 a 2 T ^3 ^ ^ '^'c 3 a 2 ~ 

(c3^a3)Ca2 


-^a3 E E '^d3d2- 

d2^ci2 {d3^a3)Ca2 

The remarkable sum rule SB 




aza2 


( 2 ) 


(3) 


0,2 03 


allows one to construct the Schrodinger equation solution T. 

DYE in configuration space become partial differential equations and re¬ 
quire asymptotic boundary conditions to define a unique solution. Let us now 
introduce the relative coordinates for the four-body problem in order to fix the 
representation of operators in Eqs. (§) and to describe the asymptotics of the 
components ^' 0302 - There exist 18 sets of relative Jacobi coordinates in one to 
one correspondence to 18 chains of partitions 0302 - Note, that among them 
the only two sets are topologically different. One corresponds to chains with 
partitions 02 of the type 3-1-1 {e.g., (123)4) and the other one to chains with 
partitions 02 of the type 2 -|- 2 {e.g., (12)(34)). From now on we will assume 
that all the particles have the same mass, say m. The relative coordinates can 


be defined trough the radius vectors of particles by formulas |]^, ^] 


Xa3 =ri- Tj 

y03,02 = (d + ri)/2 - rfc (4) 

Za 2 = (d + Tj+ rfc)/3-r. 
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for 02 = {ijk)l, as = {ij)kl and 


Xbg = Tj - Tj 

yb3,b2 =rk-ri (5) 

Z62 = + rj)/2 - (rfc + ri)/2 

for 62 = {ij){kl), 63 = {ij)kl. With coordinates (^) and @ the operators Hq 
and Va^ have the form 

^0 ~ + ^a. 2 ^ya^,a 2 /^a2 ) ’ 

= Kj 3 (Xa 3 ), 

where A stands for Laplacian and aa2 = 2 >lA, = 2/3 for the case of partitions 

02 of type 3 + 1 and aa2 = 1 , /3a2 = 1/2 for partitions 02 of type 2 + 2 . 

The asymptotic boundary conditions for solutions of DYE (^) depend on 
the problem under consideration. For the bound-state problem one should look 
for quadratically integrable solutions of Eqs. (§) satisfying 

J a^^a 2 dZa 2 I ^a 3 a 2 (^03 > y “3,02) ^<12)1 < OO) (6) 

which will give a quadratically integrable wave-function T, according to Eq. (^). 
In fact, for the short-range potentials the components ^a2,a2 have exponentially 
decreasing asymptotic behavior 0 

^Ct3^2(^) ^ -^CL30'2\^) |X|^/^ ^ 00 ). (7) 

In this equation we have introduced the shorthand notation X for the 9 -vector 
{xag, ya3,a2 > } and X for the unit vector X/|X| . For convenience of notation, 

let us define six-vectors x^j = {x^g, ya3,a2} and ya3 = {ya3,a2,202}- The vectors 
Xa2 and yag have a simple meaning. The first one describes the relative positions 
of particles belonging to subsystems of the partition 02, and the second one 
corresponds to the configuration of particles k, I and the center of mass of the 
pair (ij) for 03 = {ij)kl. Note, that for the four-body position vector X we 
have the representations 


^ — {^a3)ya3} — {^a2)^a2}- 

In the case of the scattering problem we will deal with the functions marked 
by sign +. For these functions the asymptotic form of the components 
depends on the initial state of the system. We restrict ourselves to the scattering 
with two clusters in the initial state. Two-cluster channels can be characterized 
by the bound states of the two-cluster Hamiltonians 

^ 0.2 — Ta 2 + h/3 = ~~(Axag +0^02^X03,03) T (8) 

CL3Ca2 CL3Ca2 
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describing the three-body subsystem {ijk) for 02 = {ijk)l and the two non¬ 
interacting pairs {ij) and (kl) for 02 = {ij){kl). We denote by the 

Faddeev components of the bound state of the Hamiltonian Ha 2 obeying the 
Faddeev equations 

{Ta2+Va,-ea2Wal = -ya, E > 

(c37^a3)Ca2 

where £a 2 is the respective binding energy. 

In order to characterize the three-cluster breakup channels let us introduce 
the two-body bound-state wave function ipasi^a^) being the solution of the 
Schrodinger equation 

{Has ~ ^as)''Pas — + i43(^a3) ~ ^03 }' 0 a 3 (^03 ) = 0 . 

For the sake of simplicity of notation in what follows we will assume that the 
Hamiltonians Ha 2 and Has possess at most one bound state. 

The four-body wave function components scattering with 

two clusters in the initial state can be arranged as the sum ii 

Ti|22(X,p/2) = 6a2l2(t>al{^a2) exp{ipa2 ' 

+(I^Z^^^2)Ua2l2{'^a2,Vl2) + '^as{^as)Kll2^ycis^Vh) + t^a2a3,Z2 PZ2 ) • 

The first term in (^) is due to the initial state with two clusters being subsystems 
of the partition I 2 and moving with the relative momentum p; 2 . The total 
energy E and the momentum p /2 are related by the expression 

E = 6,2 +fii2 — {pi2f. 

The second term in describes the elastic (02 = I 2 ) and the rearrangement 
(02 / h) processes. For large separations, the functions Ua 2 i 2 become spherical 
waves: 

TT t \ A exp{iv'.F-ea2|za2|} ,, I ^ 

^a2/2 (^Q'2 ; P /2 ) ^ I I ( 1 ^ 0 , 2 ! ^ (^0) 

\'^a2 \ 

with Aa 2 i 2 being the amplitudes of the processes (2 —> 2). The remaining terms 
of dH) correspond to breakup processes (2 —> 3) and (2^4). The asymptotics 
of the functions U have the form 

Jja2 ( '1 ja2 6 Xp{i\/-£i — 6^3 IYos |} ,, , _^ n H 1 i 

^ 03(2 wa3) P*2) ~ '^a3;2 |y |5/2 (ly^sl (H) 

TT ^ \ . A exp{i\/H|X|} 

k'a2as,l2\^^Pl2) ^ Aa2as,l2 |X|^ CX)j. 

Here, the amplitudes A are related to the breakup amplitudes by the formulas 

/,3,2(2-3)= E Kh 

a2Das 
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/;,(2^4)=^^A2a3,i2- 

0,2 03 


The asymptotic boundary conditions dHU) make the solution of DYE (|2|) 
unique. For the more complicated case of scattering states with three and four 
clusters in the initial state we refer to the paper [17|. 


3 Cluster reduction of the DYE 

In this section we consider the four-body bound state problem and the low- 
energy scattering problem when only the two-cluster channels are open. The 
latter means that the energy of the system obeys the inequalities 

^ - £02 >0, E - £a3 < 0, E < 0. 

Under these conditions the two last terms in (^) describing virtual breakup 
processes vanish exponentially as |X| approaches 00 and the asymptotics of 
'I'ata 2 becomes simpler. Note, that every component ^a 3 a 2 both bound and 
scattering state now is an exponentially decreasing function of the coordinates 
Xa 2 corresponding to the separation of the particles inside the subsystems of 
the partition 02 . It means, that the essential part of the component '^a 3 a 2 is 
concentrated in a tube-like domain of the configuration space such that < 
Ra 2 with sufficiently large parameter Ra 2 ■ It is natural to find the approximate 
solution of Eqs. @ by imposing the following boundary conditions: 

^^0302 |ra2 ii; 

where T^j means the boundary of the cylinder; jx^jl = Ra 2 - Imposing these 
boundary conditions is the first step of the reduction. As a second step, we 
expand the components in the following series: 

00 

^a3a2(Xa2,Za2) = Y) ■ (13) 

1=0 

Here (zaj) are new, unknown, functions and are the eigenfunctions 

of the Faddeev equations for the subsystems of the partition 02 

u* + + u, i: = (n) 

{c37^03)Ca2 


obeying the boundary conditions 


€l,li^02)\Ta 


= 0 . 


(15) 


The equations (14) with the spatially restricted conditions (15) are known to 
have a purely discrete spectrum of real eigenvalues |]^]. The functions ^ 
and the eigenvalues are assumed to be enumerated in order of growth the 
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quantities • One observes, that for large enough value of Ra2 , the functions 
(j)^^ g and the quantity approach the ground-state Faddeev components and 
the bound-state energy of the Hamiltonian respectively. The remaining 
part of basis functions extends the function (^>“2,0 ^ complete set |p^ . 

This set of eigenfunctions is not orthogonal due to the nonsymmetry of ( |l^ 
with respect to the superscript C3. The biorthogonal basis is formed |^| by the 
solutions of the adjoint to Eqs. (0): 




+ E = 4C,;.- 

(c3f^a,3)Ca2 


(16) 


€lk read 


The completeness and biorthogonality properties for the functions ^ and 

E = Kc3^(^a2 - x'aj, 


k =0 


a 3 Ca 2 ' 


Using these conditions, we express the coefficients in (H) as the projections 


FLi^a2)= E /^Xa2V’aLi(^“2)^a3a2(Xa2,Za2)- (17) 


a 3 Ca 2 


Substituting the expansions (|^ into the DYE (D and projecting onto 
the biorthogonal basis functions ^ (xa2) we get the final equations for the 
coefficients F^^{za 2 )- 

- E + = 

- E ^ ^ E'l’tjMFiM)- ( 18 ) 

a.3Ca2 d2y^a2 (rf3^a3)Ca2 ^>0 

Here the brackets (• |•) stand for the integration over Xa 2 in the domain 1x^2 1 < Ra 2 
and the coordinates Xd 2 and Zd 2 are assumed to be expressed through the set 
(^03) ya 3 ,a 2 ) ^a2 ©HI- 

The asymptotic boundary conditions for the functions (202) can be easily 
obtained from (^ and dH) by projecting according to formula ([T^ . Eor the 
bound state-problem they read 

E / dZa2\Fa2i^a2)\‘^ < OO, (19) 

fc>0“' 


and for the scattering states they have the following ‘two-body’ form 

Fa2{'^a2) ~ exp{ipa2 • } -^ I Za2 exp{i E - |Za2 |} (20) 
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for the open channels and 


k>l (21) 

for the closed channels. The equations (p^ are the desired coupled channel 
equations for four distinguishable particles. 

Let us now turn to the case of identical particles. On the first stage, note 
that due to the identity of particles, it suffices to define only two Yakubovsky 
components , namely, one for the partition 02 of the type 3 + 1 and one 
for the partition 02 of the type 2 + 2 . We fix these two components choosing 
03 = ( 12 ) 34 , 02 = ( 123)4 and 63 = ( 12 ) 34 , 62 = ( 12 ) ( 34 ) and denoting = 

Ti, = ^'2- 

Let be the cyclic and anticyclic permutation operators for four particles 
and be the same for a three-particle subsystem where the subscipt i cor¬ 
responds to the fourth particle, that is not participating in the permutation. 
With these operators the representation (^) can be rewritten in terms of the 
components Ti and '^2 as follows 

T = [/ + + P+P+ + aP-] [I + +4+ + P^]^i + [/ + P4+ + Pf ] [/ + P+P+]^2. 

The DYE (^) for the functions 'Lj, z = 1 , 2 read 

(Ho + V -E)^i + ViP+ + P^)^i = -Y[(P+ + aP+)Ti + (P+ + P4+)^2], 
{Ho + V - E)'if2 + V{P+P+)^2 = -V[P++ aP^]P+^i. 

Here V = E(xi2) and a = + 1 (— 1 ) for the boson (fermion) case. 

The expansion (p^ ) assumes the form 

00 

^i{xi,yi) = * = 1,2, (22) 

1=0 

with the basis functions cj)] being the solutions of the symmetrized Faddeev 
equations 

(Ti + vm + V{P+ + P^m = 

{T2 + V)cj)f + VP+P+(fi = e^^ 4 >l 

The biorthogonal basis functions ipf are the solutions of the respective adjoint 
equations 

(Ti + V)^Pl + (P4+ + P^)V^I;I = 


{T 2 + 1/)V’? + P^P^V^pf = 

The expansion (^) leads to the symmetrized form of (^) 



(23) 


-E + e^] F^(zi) = -{^Pi\V[P+ + aP+] ^ 

I J l >0 

-{'iPl\V[P+ + P+]Y.(^lFr), 

m>0 

j--A,, -£; + eAF|(z2) = -{'il^llViP,^ + aP+]P+J24>}Fl). 

[ rn j 

The boundary conditions for the equations (p^) follow directly from (|l^ , (p 0 |), 
and (^l|). 


4 Application to the four-nucleon system 

To apply the equations (|^ ) to the four-nucleon system one needs to perform a 
partial-wave analysis. Taking into account the conservation of the total angular 
momentum L and the total spin S of the system for chosen N — N forces S 3 
and MT I-III we can use the L — S coupling scheme and obtain for the s—wave 
components 

Fli{z) = \z\j dzFi{z), k = l ,2 ( 24 ) 

(here and in what follows the index S is introduced to distinguish the states 
with the total spin S' = 0 or S = 1 ) the following resulting s—wave equations]^ 

(“2^2 - e + £s,m)Fs,mi^) = 

/ dv - D^'^^li{x 3 ,y 3 ,)Fli{z 3 )) 

^ J-1 J-1 Xsvszs ^ 

1 pi pi ^2 

- / du dv (x4,y4)T| (Z4)), ( 25 ) 

2 J-1 J-1 XiyiZA “ 

(“2^^ “ £ + e|,n)-^l,n(^) = 

[ dv Ci J2^h(^5,y5)Fliiz5)). 

J-1 x^y^z^ “ 

Here the brackets (-I-) mean the integration over x and y in the domain x‘^+y‘^ < 
E? and summation over the components of vector-functions ^ and ^ being 
the solutions of s—wave Faddeev equations 


^ We have restricted the summation in formula (^) to hnite numbers Ni and A2, so that ( 25 ) 
should be considered as approximate equations,which become exact in the limit Ai, N2 —> 00. 


9 







{-dl - \dl + Vf{x)\f dv^^Bf<^>li{xi,yi) = 
14" J J-i xiui 

{_a2 _ dl + $|^„(x,y) + F/(x)i?fc^|,„(x2,y2) = (26) 

= 4,n^S,n(a;>y)> 

and adjoint equations 

|-4 - 7^^ + + / dv^^Bfvf {xi)'i/l^{xi,yi) = 

I 4 " J ’ J_i xiyi 

= 4,z^4(^’y)’ 

{_02 _ q2 ^ + B^Viix2)^l^ix2,y2) = (27) 

= 4,n^l,n(a^>y)- 

In the above equations we have renormalized the energy of the system according 
to e = mE/h?. The matrices B^, and C^, /c = 1, 2 realize the representation 
of operators , P+P+, P^ — P+, P+ + P+, and (P+ — P+)P+ in the 

spin-isospin space of the four-nucleon system, respectively. We give their values 
for the systems nnpp and nnnp in the Appendix. There we have described the 
representations of interaction potentials 14 ^ and values of coordinates Xi,yi,Zi 
appearing in the equations. 

The components Fgi{z), due to (^^, obey the regularity conditions 

FiM = 0 . 

The asymptotic boundary conditions for the bound-state problem according to 
(HD and (p^ have the form 

poo 

Y, dz|P|/z)|2<oo, A: = 1,2. (28) 

z > 0-^0 


For the n-'^H scattering problem the only three-body subsystems have the 


bound state, so that (|^) and {21) lead to the conditions 


Pgg(z,p)~p ^ sinpz-|-a 5 (p) cospz, 


Fli{z,p)^Q, l>l, (29) 

Fs,iiz,p) ~ 0 , I >0. 


For the zero relative n-^H scattering energy the first condition of (29) should 
be replaced by 

fUz,0)^z-As, (30) 
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where As is the n-^H scattering length. The relationship between the scattering 
amplitude as{p) and the phase shift 5s{p) is of the form 


as{p) = 


tan (5s (p) 


P 


(31) 


Alternatively, the scattering length As can be obtained from as{p) by taking 
the limit, 

As = - linios(p). (32) 

p —>0 

The equations ( pSD were solved numerically by making use of a finite- 
difference approximation for the differential operator and a spline expansion 
for the integrands on the right hand-side. The mesh parameters were varied 
extensively to obtain the relative numerical uncertainty less than 1%. The ba¬ 
sis functions and ^si were constructed numerically as the solutions of the 
eigenvalue problems (p^ and (^). The maximal nonorthogonality of the basis 
functions caused by the numerical solutions of 


and (27) was given by 


We have used S3 and MT I-III s—wave potentials for the bound-state cal¬ 
culations. For n-^H scattering we have used only the MT I-III potential. These 
potentials are known to have a strong infinite (MT I-III) and finite (S3) repul¬ 
sive core. Such a singular behavior of the potentials requires an appropriate 
treatment of the problem, so we have chosen these potentials to demonstrate 
the efficiency of our approach. On the other hand, the data of the direct (ex¬ 
cept for numerical approximations) solutions of the four-nucleon Yakubovsky 
equations (YE) for bound states with these potentials are available |^, IT | which 
facilitates the comparison of the results of calculations. 

In Table 1 we show the convergence of our calculations for '^He binding 
energy with respect to the number of terms taken into account in the expansion 


(22) and in Eqs. (25), respectively. As one can see a converged result within 
the uncertainty less than 1% can be reached with about 10 basis functions. In 
Table 2 we have collected our final values for ^He binding energy and the data 
of solutions of YE gsni. The agreement of calculations is near perfect. 

The calculations of n-^H scattering were performed for the zero neutron 
laboratory energy. We have used the boundary conditions (^) for zero energy 
calculations. The scattering length As was extracted from the asymptotics 
of the solutions according to the formula (^). We have collected the results 
of our calculations for spin-singlet and spin-triplet scattering lengths Aq and 
Ai together with available data of the other authors and experimental values 
in Table 3. Comparison of the results is hampered by the fact that different 
approaches and different models for N - N forces were used. Nevertheless, the 
agreement of the results is reasonable except for the case of the paper |23] where 
the Resonating Group Method (RGM) has been used. In our notations, the 


RGM ansatz corresponds to taking into account only one basis state in (22). 


II 









We have performed the calculations within this condition and obtained 3.41 fm 
for the spin-singlet scattering length which is very close to the result of |23|. 
Hence, we can suppose that the difference between our converged result and 
the RGM one is due to the approximations of RGM. 


5 Summary 


We have described the cluster reduction method to treat the four-body bound- 
state and scattering-state problem. The method leads to GRG equations for 
the functions depending on the relative coordinates between the subsystems of 
the two-cluster partitions. After a suitable partial-wave decomposition, these 
equations become the one-dimensional integro-differential equations which can 
be solved numerically on rather small computers in contrast to the original 
DYE, for which numerical solution even for simple N — N forces requires su¬ 
percomputer facilities. The results of calculations given in the previous section 

for the numerical solution of 


show the efficiency of reduced DYE (23) and (E 
both the bound-state and scattering-state problems for the four-nucleon sys¬ 
tem. The method proposed can be extended both to the case of more realistic 
nucleon-nucleon potentials in four-nucleon system, and to systems consisting of 
more than four particles. In the latter case the N-body differential Yakubovsky 
equations should be taken as a starting point for the reduction procedure. 

The efficiency of the cluster-reduction method for the scattering problem 
demonstrated for the n-^H system suggests it will be useful for the direct four- 
body calculations of various multichannel reactions. 
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Appendix 

In this Appendix we describe the representation of the potentials and the 
matrices , and C^, k = 1,2 and give the formulas for the coordinates 

Xi,yi,Zi appearing in the integrals of Eqs. (pB]), (|2^), and (p^). For the bound 
state of the nnpp system the total spin equals to zero, 5 = 0. The operators 
are the diagonal matrices such that 

= V^ = diag{v\v^}, 

where and u* stand for the singlet and triplet N — N potentials, respectively. 
The matrices B^, and have the form 


= D^ = C? = 



dO _ ^0 _ 

^2 — ^2 — 

For the system n-^H the values of the total spin 5 equal 0 or 1. The 
operators are the diagonal matrices 

if = diag{v^,v\v^}, 

V,^ = v^, 

Vi = diag{v^ ,v^ ,v^ ,v^}, 

V2 = diag{v^ . 

The explicit form of the matrices B^, , and is as follows 


50 = 


/ i 

“4 

V 0 


( 

1 

4 

0 -f 


B^ = -1 


Bl = 


1 

4 

3 

4 
0 


= 


0 


0 \ 
0 
0 


Vo 0 0 -IJ 
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= 


-T2 i 

1 - 1^2 
i \/2 0 -i 0 


V 0 1^2 


0 


-u 


cV = 


(IV2 

kVQ 

V 


Cl = i-^V2 0 -l/Vs), 

0 


c} = 


( -j^^/3 i^/3 

1^/3 


V 


4 V - 12 

0 -i \/6 


_ / 0 -i^/3 0 i^/6 

- I -1^3 0 i ^/6 0 


3 V vv 3 

The expressions for the coordinates Xi,yi, Zi can be given by formulas 

\ 1/2 / Q „ 1 „ .t \ 1/2 


xi = ( - xyv 


/9 2 1 2 3 V 

) ’ = V le"" + 4^ + 


X 2 = y, y2=X, X3 = X1, X4 = Xi, X5 = y, 


n 2 , 2 

2 \i/2 

1 , 

2 , 2 

16 \ 

y 3 = f gyi + 

1 

^3 = 3 ! 

- 

yyi^nj 


^4 2 2 4 \T /4 2 1 2 2 \ 

y4 = ( gyi + + -^yizuj , ^4 = (-yi + -z - -yi2;u 1 

ys = Qx^ + 2:^ - xzx^ , 25 = ^ (^x^ + 2^ + 2 xzv^ 


1/2 

\ 1/2 
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N 

S3 

MT l-lll^ 

MT 1-1II2 

1 

23.86 

24.6945 

25.154 

2 

27.11 

28.3060 

28.869 

3 

27.23 

28.4734 

28.975 

4 

27.75 

28.7086 

29.223 

5 

28.38 

29.3345 

29.688 

6 

28.42 

29.4204 

29.783 

7 

28.39 

29.9648 

30.284 

8 

28.55 

29.5980 

30.277 

9 

28.67 

30.0534 

30.304 

10 

28.72 

30.0829 

30.315 

11 

28.71 

30.0901 

30.307 

12 

28.74 

30.1369 


13 

28.73 

30.1392 



Table 1. The convergence of the ^He ground-state energy i5^(MeV) with 
respect to the number = A^i = Ai 2 of terms in (^) taken into account in the 


numerical solution of Eqs. 
MT I-III^ from g 0. 


Here MT I-III^ is the potential from 0 and 


17 





















Refs. 

S3 

MT I-III^ 

MT 1-111 2 


i 

28.7843 


30.3117 

1 


28.80 


30.29 

1 

g 


29.6(2) 


present work 

28.7(4) 

30.1(3) 

30.3(1) 


Table 2. The results of the calculations of the ^He ground-state energy 
£^(MeV) for spin-dependent S3 and MT I-III N — N forces. Here MT I-III^ is 
the potential from |Q and MT I-III^ from [P, pT| . 
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Refs. 

Aq fm 

Ai fm 

Present work 

4.0 

3.6 

0 

4.09 

3.61 

0 

4.23 

3.46 

H 

3.905 

3.597 

n 

3.8 

4.9 

n 

3.38 

3.25 

n 

4.13 

3.73 

exp. |24] 

3.91 ± 0.12 

3.6 ± 0.1 


Table 3. The calculated spin-singlet and spin-triplet Ai scattering 
lengths using the MT I-III^ potential of |]lj] , the results of other authors and 
experimental data for n-^H scattering. 
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